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Abstract. In this paper the angular characteristics of a multimode fiber SPR sensor 
are theoretically investigated. By separating the contributions of beams incident at 
different angles, a compact model is presented to predict the shift of the resonance 
wavelength with respect to the angle and the environmental refractive index. The result 
suggests that the performance of conventional fiber SPR sensors can be substantially 
improved by optimizing the incident angle. Furthermore, our investigation suggests 
some problems in previous reports. 
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Recent applications of fiber SPR (Surface Plasmon Resonance) sensors to biochemical 
analysis have attracted much interest m ei m a un m el m- By integrating SPR 
region on the cylindrical surface, a compact multimode fiber SPR sensor is produced. 
Fiber SPR sensor works based on attenuated total reflection (ATR) spectroscopy. For 
wavelength interrogation, a modification of surface mass load on the fiber sensor will 
change the resonance condition, which results in a spectral shift of the resonance dip. 
Compared with other variants, multimode fiber SPR sensors are simple, robust and 
easy to produce mmm- To provide for a growing need in this developing held, we 
built a sensor system and created software for real-time measurement m In such a 
conventional fiber SPR sensor, broadband light in various incident angles result in an 
asymmetric absorption curve, which requires high order polynomial fitting to resolve 
the resonance wavelength.This implies that it is possible to improve sensor performance 
by restricting the angular distribution of input light to a narrow range. 

Incident angle is an important parameter in surface plasmon resonance. It 
determines the projection of the wave vector at the interface. A multimode fiber with 
a diameter of hundreds of microns is rigid enough to keep straight and to preserve the 
propagation angle of a light beam. This is the experimental basis of our investigation. 
Although there are already many theoretical and experimental studies of multimode 
fibers m mi ng, an analysis to distinguish the contributions of light beams with 
various incident angles is still missing. Meanwhile, the community lacks a clear picture 
for sensor prediction and optimization. This motivates us to re-visit this subject from 
a theoretical perspective. In this paper we focus on the angular characteristics of the 
multimode fiber SPR sensor under wavelength interrogation. The main framework is 
based on Kreschmann’s theory and Fresnel reflections in multilayer films. 

2. The calculation model 

The schematic of a multimode fiber SPR sensor is shown in Fig. [l] The light is launched 
through a ring mask for angular control. Collimated beams are then focused at the center 
of the input face of the fiber with a microscope objective. Assuming the fiber sensor is 
fabricated from a plastic cladding silica fiber with diameter d = 600 p/m, the removed 
cladding length is L = 4 mm and the thickness of the gold him is t = 45 nm. 

Considering a circularly polarized light source, its transmission as a function of 
incidence angle and wavelength is given by mm- 



(i) 


where N = L/d tan 6 is the number of reflections. The transmitted light is divided into 
two components, with p- and s-polarization. The quantum factor of the CCD is not 
included since beams are considered in one angle with same efficiency. 
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Figure 1. A multimode fiber SPR sensor with controlled light angle. 


The cascaded reflectance in three-layer films is given as follows [13] : 


R = 


r me exp (2 ik mz t) 


( 2 ) 
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where the fiber-core, metal and environmental layers are denoted by c, m, and e, 
respectively. The fiber-core/metal and metal/environment interfaces are labeled as ‘cm’ 
and ‘me’. The reflection coefficients with s- and p- polarized light on the interface are 
given for wave vectors along the z-axis and dielectric constants mm- 
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The refractive indices of the fiber core are interpolated from [To], while the dielectric 
constants of gold are according to Johnson and Christy HU. Wave vectors in different 
mediums can be written as: 

h z = (£i~ - k 2 x ) 1/2 (7) 


k x = Hcrr-sinO 
An 


( 8 ) 


where c is the speed of light and Ao is the wavelength of the incident beam in free space. 
The letter i can be c, m, or e depending on the medium in question. k x is the component 
of the wave vector parallel to the interface, which is conserved in all media. In principle, 
the incident angle of the light beam 6 is bounded by the critical angle 9 cr = sin ~ 1 (n c i/n c ) 
and 7t/2. The fiber numerical aperture NA = \Jn 2 c ~ rfy is also related to these indices. 
For a multimode fiber with a numerical aperture NA=0.50, the critical angle is 69.8° if 
n c = 1.446 RIU is used. In our investigation, we discuss incident angles ranging from 
65° to 90°. 
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Assuming the fiber SPR sensor is immersed in the water (n e = 1.33 RIU), Fig. [2] (a) 
presents a three-dimensional plot of the transmission ratio T = T(1.33, 9, A) with respect 
to both the angle and the wavelength. For specific values of the incident angle (9 ranging 
from 72° to 87° with an interval of 3°, as well as 89°), black lines outline their absorption 
spectra. The absorption peak positions (corresponding to minimum transmission) are 
traced by a cyan line to show the overall tendency. Clearly, decreasing the incident angle 
increases the resonance wavelength. The reason for this trend is not straightforward, 
since the surface plasmon vector k sp changes as the light wavelength increases. To 
explain this, we define an auxiliary parameter /3(Ao) = £ m (Ao)£ e /(£m(Ao) + £ e ), which is 
a function of the wavelength Ao- To satisfy the resonance requirement, the component 
of the wave vector along the interface should be equal to that of the surface plasmon: 


I3(\ 


oj 


(9) 

( 10 ) 
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Thus the resonance position is defined as a relationship between the angle and the 
wavelength. 


9 = sin 1 (y'/3(A 0 )/n c ) (11) 

In Fig. [2] (b), the blue line presents resonance points obtained from Eq. 11, while the 
cyan line is the locus of peak absorption from Fig. |2](a). The two curves do not perfectly 
overlap, since they arise from different models; the blue line is derived assuming an 
ideal interface between two semi-infinite media, while the cyan line represents a working 
sensor with a thin film. The discrepancy between the two models can be accounted 
for by invoking an ‘effective’ dielectric constant. As the metal film thickness increases, 
the negative effective complex dielectric constant £ e g on the metal side of the interface 
decreases, which is equivalent to the increment of n e for they are equal in the expression 
of {3. Thus the resonance wavelength increases slightly and the cyan line moves toward 
the blue line. In theory, the sensor interface approaches to the ideal model as the film 
thickness increases. 

To investigate the performance of the fiber sensor in different environmental 
conditions, resonance lines for n e from 1.33 RIU to 1.40 RIU are stacked in Fig. [2] 
(c). The z-axis is the environmental refractive index. The lines are colored according to 
the transmission ratios at the resonance positions, as indicated by the color-bar. The 
heavy magenta lines are contours of constant incidence angle. These same magenta 
lines are plotted in Fig. [ 2 ] (e) for clarity. I11 Fig. [2](f ) we present the corresponding 
sensitivities, which is defined as the ratio of the shift of the resonance wavelength to the 
change of environmental refractive index. For our case, we fit the discrete series first, 
and then calculate the corresponding value of the derivative of the fitted function. The 
meaning of the heavy black line is discussed below. Furthermore, Fig. [2](d) presents the 
transmission ratios for two specific values of the incident angle, 9 — 81 and 9 = 75, to 
give a direct impression of their shapes and FWHMs. 
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After processing, we achieve a systematic understanding for the performance of a 
multimode fiber SPR sensor. In particular, Fig. |2](c) presents an elegant perspective for 
sensor optimization, and from it we can draw three conclusions. First, the working angle 
is limited; it ranges from 69° to 90° in our setup. Therefore, there is no justification 
for using a base fiber with a numerical aperture larger than 0.5. Second, lower angles 9 
give higher sensitivity. The tradeoff is that its measurement range is narrower than that 
with large angles. To obtain high sensitivity, we should move the working position to a 
magenta line towards the right side of Fig [2] (c). For example, when n e = 1.35 RIU, the 
sensitivity is 2144 nm/RIU for a light beam with angle 9 = 84° , while for a beam 
with 9 = 75° (near the edge of a multimode fiber with NA=0.39), the sensitivity 
can be improved to 7101 nm/RIU, almost triple that for 9 = 84°. These values 
are highlighted as red points in Fig. [2] (f). For an even higher environmental index 
n e = 1.37 RIU, the corresponding sensitivities are 17134 nm/RIU and 3719 nm/RIU, 
differing by a factor of five. Third, when the selection of the incident angle is adapted 
to the environmental refractive index range to be measured, the maximum sensitivities 
under different conditions are close to each other since their working positions are located 
on the right edge of the spherical face. For n e = 1.39 RIU, the corresponding sensitivity 
is 17832 nm/RIU if a light beam with 9 = 78° is used. This result is close to the red 
point of 17134 nm/RIU under n e = 1.37 RIU in Fig. [2] (f). 

To compare our analysis with previous investigations, we model a multimode 
fiber SPR sensor with incoming light comprising a distribution of incident angles, 
corresponding to a numerical aperture NA=0.22. This model is henceforth referred 
to as a ‘conventional fiber sensor’. Other parameters are the same as for the previously 
discussed single-angle model. Transmission ratios and resonance wavelengths for 
different n e are calculated, taking into account the angular dependence of the power 
distribution and the CCD efficiency [13, 0 T3]. To allow direct comparison, we 
introduce a concept of “effective angle”, whose value is interpolated from the resonance 
wavelength of the corresponding resonance line. Thus we project the result of the 
conventional fiber sensor onto the spherical face in Fig. [2] (c), as illustrated by the heavy 
black line. As shown, the black line lies almost on top of the magenta line corresponding 
to 9 = 84°, so the effective angle of the conventional fiber sensor is about 9 = 84°. The 
black line lies entirely in the blue (high absorption) region, so the red (low absorption) 
zone near the edge of the plot, which corresponds to light beams in angles greater 
than 87°, is inactive in practice. It suggests that light beams being almost parallel to 
the fiber axis, corresponding to the red edge of the surface in Fig. [2] (c), contribute 
little to the sensor performance because of their weak absorption. As n e increases, the 
resonance wavelength increases and the wave vector k x = k sin 9 decreases; this leads to 
preferential weighting of smaller angles, which should account for the deviation of the 
black line towards a smaller effective angle at large refractive indices in Fig. |2](c). 

The transmission curves for a conventional fiber sensor for environmental indices 
n e = 1.33 RIU, n e = 1.365 RIU, and n e = 1.40 RIU are plotted in Fig. [2] (d) with legend 
‘NA’. The resonance absorptions are less than 0.4, which suggests a moderate resonance 
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state caused by the short SPR region. In general, transmission for the conventional 
fiber sensor is a superposition of resonance curves of various angles. For the condition 
n e = 1.40 RIU, the FWHM of the spectrum is close to 200 nm and the bottom of the 
curve is nearly flat, which makes it hard to identify the minimum position. Compared 
to the red lines (6 = 81°), the FOM of the sensor is much poor. In Fig. [2] (e) and (f) 
we show resonance wavelengths and sensitivities of the conventional fiber sensor (black 
line and star markers). These results are somewhat higher than for light with a single 
incidence angle of 6 = 84°. The sensitivity is 1523 nm/RIU at n e = 1.34 RIU, as labeled 
in (f). This value is close to what we found in our experimental study [T2j, if we neglect 
some differences in sensor parameters. 

Except the low efficiency and moderated performance in the conventional 
configuration, many earlier reports did not take a serious consideration for the angular 
characteristics of the fiber SPR sensor m m ei For example, the Y-combiner, 
which is widely found in the literature and involves thermal tapering and then fusing 
together of fibers, will change the angular distribution of the propagating light and the 
re-distribution in angle is out of control. Therefore it is impossible to trace the paths 
of rays. But SPR effect strongly depends on the incident angle. Thus it leads to a 
mess in theory. More attention should be paid to this issue. Based on our investigation, 
we recommend that researchers working on multimode fiber SPR sensors stick to the 
original configuration without Y-combiner, although it is a bit inconvenience. Another 
important problem is the mismatch in aperture size between the fiber itself and the 
spectrometer, which had been neglected in both the experiment and its theoretical 
simulation. For compactness and convenience, most fiber SPR configurations use a fiber 
spectrometer to acquire spectra. However, the fiber spectrometer has its own numerical 
aperture, typically NA = 0.22. Obviously, under this configuration the large-angle 
contributions are wasted, even though the base fiber has a large NA. A possible solution 
is to use another microscopic objective to recover a quasi-collimated light beam before it 
enters the spectrometer; otherwise one could use a customized fiber spectrometer with 
a larger NA. 

About twenty years ago several fiber SPR sensors with oblique light input were built 
for angular interrogation [HHIDS) 2D]- Obviously, their configurations were much different 
from our proposal as we use white light and wavelength interrogation while laser and 
angle shift are adopted in those reports. The kernel of our proposal is to purify the input 
and evaluate the contribution separately. This can not be achieved in these mentioned 
papers where non-meridian light beam dominate the effect. In addition, if the ring mask 
is thick enough, the skew rays should be reduced greatly in our configuration. 

The number of reflections N ( 6 ) is determined by the fiber diameter and the sensing 
length. Since power operation doesn’t change the position of the minima, these two 
parameters have little effect on our model if the absorption threshold hasn’t been 
reached. This is also demonstrated in Fig. [2] (b). The red line, which represents the 
result of a single reflection of R p in Eq. ([!]), overlaps almost exactly with the cyan line 
of minimum transmission; the small difference is a numerical artifact arising from the 


Angular characteristics of MM-fiber SPR sensor 


7 


minimum search algorithm. In fact, if the averaging reflection times is fixed as N=l, the 
SPR model reduces to a prism SPR model. In this case Fig. [2] (c) maintains its shape, 
and corresponding conclusions still hold. Thus our model may also serve as a reference 
for prism SPR studies. 

When the sensing length increases, N ( 9 ) increases, the resonance is intensified and 
the red zone in Fig. [2j^c) shrinks to the edge. Consequently, the importance of light 
beam with large angle increases; working positions of the conventional fiber SPR sensor 
move toward the edge and the sensitivity decreases. This should account partly for 
differences between our numerical results and previous reports HHng. 

As illustrated in Fig. |2][b) , the him thickness make a slight difference between a fiber 
SPR sensor of finite thickness and the ideal model. If the him thickness changes, the 
cambered surface in Fig. [2] (c) shifts slightly. Furthermore, the sensitivity is determined 
as a slope of the magenta curves. Thus the variation in performance is delicate and 
it is difficult to identify a general tendency. In contrary, the him thickness plays an 
important role in the signal amplitude, which should be considered before other factors. 

Fiber SPR tips or tapered hber SPR can also be described qualitatively in our 
model. As proposed in works [21, [22], post-processing of the hber will directly 
decrease the incident angle of the main portion of the light with respect to the normal 
of the cone surface. Thus the effective working angle decreases and the sensitivity 
is improved. Meanwhile, our model helps to identify the maximum angle for post¬ 
processing: according to Fig. [2] (b), the tapering angle should be less than 20° under 
environmental index n e = 1.33 RIU. 

4. Conclusions 

In conclusion, we have numerically calculated and analyzed the angular characteristics 
of a multimode hber SPR sensor. A general model is presented for sensor prediction and 
optimization. It suggests that small-angle components of the light beam have better 
sensitivity. We have proposed an angular hber SPR sensor based on a multimode hber. 
For a specihc measuring index, the performance of a multimode hber SPR sensor can 
be substantially improved by deliberate control of the beam angle. 
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Figure 2. (a) Three-dimensional plot of the transmission ratio as a function of both 
wavelength and incident angle, for n e = 1.33 RIU. (b) The ideal resonance requirement 
(blue), the general trend from the full model (cyan), and the result for a single 
reflectance R p (red), (c) The stacked resonance positions for n e ranging from 1.33 RIU 
to 1.40 RIU. The black dots are the projected working position of a conventional fiber 
sensor, (d) Sample spectra for beams with incident angles 9 = 81° (red) and 6 = 75° 
(blue), as well as the conventional fiber sensor (black), (e) Resonance wavelength 
versus environmental index for a range of specific values for the incident angle, (f) 
Sensitivity versus environmental index for the same incident angles as (e). 






























































